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ABSTRACT 



We present a study of the ionization and thermal structure of neutral hydrogen clouds located 
in the Galactic halo, immersed in the extragalactic background radiation field, and supposed to 
be in pressure equilibrium with the surrounding (presumably) hot medium. The problem is solved 
numerically, but an useful analytical approximation has been derived as well. We discuss the two 
main parameters of the problem, i.e., the radiation field and the halo pressure, in the framework of 
the current models for the halo/disk interaction. We find a well defined relation between a critical 
column density at which the cloud starts to develop a cold central core with the cloud linear size. 
Making use of this relation, we suggest a straightforward method to derive the distance to the cloud. 
We discuss the possible sources of error in this determination and it is found that the method is 
particularly suitable for those clouds which are subsonic with respect to the surrounding medium, 
while for the supersonic ones, the method can only give a lower limit to the distance. The Ha 
emission from the partially ionized edge of the cloud is calculated and compared with the available 
observations; this measure is a powerful indicator of the possible presence of a shock in the cloud. 



1. INTRODUCTION 



One of the most long standing problems in the study of the interstellar medium is represented 
by the so-called high velocity clouds (HVCs). Generally speaking, HVCs were hysterically defined 
as clouds of neutral gas with a velocity relative to the LSR |fL5/?| ~ 90 km s _1 . Here we use 
the term in a somewhat looser manner, since we mean the HI cloud with velocities incompatible 
with a simple model of galactic rotation and structure. The only requirement is that the clouds 
are located well above the main Galactic disk into the halo. The mass flux on the galactic disk 
provided by the HVCs and the VHVC, very high velocity clouds, is about 0.5M Q yr _1 (Mirabel 
1989); Sembach et al. (1991) and Wakker (1991) derive a value of about 5M Q yr _1 including all 
HI with velocities not fitting galactic rotation. An estimate of their metallicity given by de Boer 
& Savage (1984) is Z < 0.3Z Q ; similar results are obtained by Blades et al. (1988) for two HVCs 
in the direction of SN1987A. HVCs are usually observed in HI 21 cm line emission; however in 
the last years a number of different observations have been attempted. Colgan et al. (1990) (and 
references therein) have reported a negative result in detecting the 21 cm absorption line, which 
they interpret as a requirement for an external heating source. Wakker &: Boulanger (1986) have 
used the IRAS data to search for 100/um emission from dust inside the HVCs but none of the sample 
clouds has been detected. This implies a lower-than- normal dust content or a low temperature of 
the dust. There are no direct evidences of molecular hydrogen related to the HVC, even if Rohlfs 
et al. (1989) suggest that the high latitude molecular cloud G90+38 can be spatially connected 
with an infalling HVC. Absorption lines towards halo stars and extragalactic objects may provide 
useful information on the composition of the HVCs. However the only firm detection of a HVC in a 
stellar spectrum is towards the star HD 135485 (Albert et al. 1989). Very recently, Bowen & Blades 
(1993) have detected for the first time Mg II absorption lines from two HVCs, which are part of the 
Complex C, in the line of sight of Mrk 205. Their results suggest a very different column density 
ratio N(Mg II) /N(HI) in the two clouds, which might be explained by a real difference in the metal 
content of the two clouds or by a difference in the ionization conditions: it is important to note at 
this regard, that the HI column densities for the two clouds differ roughly by a factor 10 (typical 
HI column densities for HVCs are in the range 10 18 < Njj < 10 21 cm -2 ). A relevant feature, 
discovered through the high resolution synthesis maps (Wakker & Schwarz 1991), is the presence of 
bright, dense cores with Nh well in excess of 10 20 cm -2 . A type of structure for the clouds made 
by bright concentrations with low velocity widths (~ 7 km s _1 ) and extended envelopes with much 
larger widths (~ 23 km s _1 ) was already indicated by Cram & Giovanelli (1976). 

As can be realized from the previous discussion, our understanding of the HVC phenomenon is 
quite fragmentary and incomplete. Definitely, the major source of ignorance is the lack of a precise 
distance estimate. The central importance of the distance knowledge has been stressed by many 
authors (Van Woerden 1993 and references therein). From the observational point of view, the 
claim made by Danly, Albert & Kuntz (1993) of the bracketing of the distance to a HVC in the 
Complex M in the range 1.5 < z < 4.4 kpc, provides at least the first firm evidence that the HVCs 
are located in the Galactic halo. The distance estimate is particularly important because most 
of the physical parameters of the clouds scale with the distance A: the physical size, £ cx A, the 
density, n oc A -1 and the mass, M oc A 2 . In addition the related problem of the origin of the HVCs 
could also strongly benefit from such an information. The most common explanation for the origin 
of the HVCs is the so-called galactic fountain model (GF) worked out by Shapiro & Field (1976) and 
subsequently detailed by Bregman (1980); Houck & Bregman (1990) extended it to the low regions 
of the halo; Li & Ikeuchi (1992) have investigated their formation in giant halos. In this model the 
HVCs correspond to the condensation mode of a thermal instability occurring in a flow of hot gas, 
generated by supernova explosions, rising from the disk. Altough the GF can explain most of the 
data concerning the HVCs (but not all: see the discussion in Wakker 1989), many doubts are still 



present on the reality of the entire circulation process. In fact, it is not completely clear that the 
superbubbles are able to breakout of the disk when the thick ("Lockman") exponential layer of the 
HI disk is considered (Mac Low & Mc Cray 1988; Mac Low et al. 1989; Tenorio-Tagle et al. 1990), 
and particularly when the effects of a magnetic field inhibiting the growth of the supernova shock- 
sweeped region are taken into account (Tomisaka 1990; Ferriere et al. 1991; Shapiro & Benjiamin 
1993; Norman 1993). On a different basis Ferrara & Einaudi (1992) pointed out that, under the 
regime prevailing in a fountain flow, dynamical instabilities, leading to convective motions rather 
than non-gravitational condensations, may have a faster growth rate, thus quenching the cloud 
formation process. On the other hand, clear evidence of hot gas (7\ ^ 10 6 K) located in the 
Galactic halo comes from the ROSAT shadowing experiments towards the Draco cloud (Burrows 
& Mendenhall 1991; Snowden et al. 1991; Herbstmeier et al. 1993). Also, Herbmeister et al. 1993 
report an enhancement of the soft X-ray emission near some HVCs in the northern sky. This result 
brings fresh support to the existence of an extended hot halo and poses some constraints on its 
physical characteristics. It is probably redundant to stress at this point the value of the distance 
determination in order to make some progress on the origin of the HVCs: if they are located in the 
halo they must be in pressure equilibrium with the hot gas in order not to be rapidly dispersed; 
furthermore the study of this interaction may lead to a better comprehension of many different 
phenomena (accretion of material onto the galactic disk, halo and disk star formation, Galactic 
chemical evolution, general gas circulation, structure of the halo) that are still unclear. 

Finally, if HVCs are located well above the main gaseous Galactic disk, they are exposed to 
the ionizing extragalactic background radiation field (EBR). For the purpose of obtaining a deter- 
mination of this field at redshift z = 0, the HVCs are an ideal "test particle", being a relatively 
quiet and little contaminated environment. This point have already been stated by other authors 
(Cowie & McKee 1976; Cowie & Songaila 1986; Songaila, Cowie & Weaver 1988); in addition, the 
local value of the EBR has been recognized to be crucial in constraining many cosmological models. 
The principal response of the clouds to the ionizing radiation is the Ha recombination emission. 
A number of searches of this emission have been performed in the last years: at least one positive 
detection has been reported (Kutyrev & Reynolds 1989). 

In this paper we are proposing a new method to derive the distance to the HVCs based on a 
photoionization model for the clouds, supposed to be irradiated by the EBR. We will show that the 
observed core/envelope structure of many HVCs can be used to derive the cloud linear size, and, 
from the knowledge of the angular size, the distance to the cloud can, in principle, be obtained. 
Of course, there are posible sources of error in the method either of theoretical and observational 
nature, and we discuss the first ones and point out the possible observational difficulties. Also, 
matching the Ha emission deduced from the model to the observed one, further insight on the 
physical state of the cloud can be obtained. 

The structure of the paper is as follows. In § 2 we describe the adopted photoionization model 
and discuss its relevant parameters; § 3 is dedicated to the results and to the distance estimation 
method, while in § 4 some additional implications are discussed. 



The equations governing the steady-state ionization and thermal balance of a cloud immersed 
in an isotropic radiation field of intensity I v are: 



n(A l ) / -Z-a^X^du [l + 0(X 1 )} + lc {X\ T)n(X>(e) = n(X i+1 )n(e)a(X i , T), (2.1) 



2. IONIZATION MODEL 
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These equations express the ionization and thermal balance, and charge conservation, while eq.(4) 
is the equation of state for the gas; X denotes the element considered and i its state of ionization; 
J u is the first moment of the field, vl x indicates the ionization limit for each species, a is the 
photoionization cross-section, <j> is the secondary ionization rate, 7 C is the collisional ionization 
coefficient, C indicates the appropriate volume cooling rate, and a is the total recombination 
coefficient. We have adopted the "on the spot" approximation in which the diffuse field photons 
are supposed to be absorbed close to the point where they have been generated. 

The various elements are divided in "primary", which enter the ionization-thermal equilibrium 
equations (2.1)-(2.4) and "secondary" which do not. In the following, H and He are considered 
as primary, while other elements (C,N,0,Si,Fe) just contribute to the cooling function and to the 
electron density (C only). Secondary elements are considered to be completely ionized by the UV 
field below 13.6 eV; a metallicity Z = 0.25Z Q and a helium abundance equal to 0.1H as been 
assumed throughout the paper. Double ionization of He has been neglected, and He fractional 
ionization has been supposed equal to the H one, x; the ionization cross section for He has been 
taken from Brown (1971). Helium is not of special importance in the temperature range that we 
consider; this justifies the rough approximation adopted. Secondary ionization rates and fractional 
heating for H and He have been taken from Shull & Van Steenberg (1985). The following processes 
have been included in the calculation of the cooling function: i) free-free from all ions; ii) H and He 
recombination; iii) electron impact ionization of H and He; iv) electron impact excitation of H and 
He (n=2,3,4 triplets); v) He dielectronic recombination; vi) electron and H impact excitation of 
secondary elements; excitation of the metastable levels by electron impact are also included. The 
obtained cooling function is almost identical to the one given by Dalgarno & McCray (1972). In 
additon to photoionization, heating is also provided by C ionization (2 eV per ionization), assuming 
that the recombination coefficient a for C is equal to the H one (Spitzer 1978). The numerical values 
for the various coefficients have been taken from Black (1981) and Dalgarno & McCray (1972). 

In order to obtain the spatial ionization and thermal structure of the cloud, eqs. (2.1)-(2.4) must 
be solved simultaneously with the radiation transfer equation. We have postulated that the cloud 
can be modeled as a slab of gas of thickness I illuminated on both sides by the (isotropic) radiation 
field I„. With the above assumption the transport equation reads 



where E 2 is the exponential integral function, I„ t0 is the field intensity at the cloud edge, and 
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which gives the solution 
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The system of equations (2.1)-(2.5) can be numerically solved trough an iterative scheme. We 
have requested the pressure P to be constant throughout the cloud; this condition is particularly 
suitable to describe the HVCs, which are believed to be in pressure equilibrium with a hot external 
medium. This argument, initially introduced by Spitzer (1956), is based on the fact that the 
pressure of the disk intercloud gas at 1 kpc would be far too low to confine the cloud, which 
would rapidly expand decreasing the density to values not allowing the detection of any spectral 
absorption feature toward distant stars. 

In analogy with the usual HII regions, we expect to find a critical length, £ c , analogous to the 
usual Stromgren radius. The latter can be dimensionally evaluated equating the number of available 
ionizing photons to the number of recombinations inside a sphere, but in principle one can calculate 
the structure of the partially ionized zone at the interface between the cloud and the surrounding 
medium using some analytical approximations for eqs. (2.1)-(2.5). If we neglect secondary and 
collisional ionizations, the electron contribution due to secondary elements and helium, one can 
write the ionization equation in the isobaric case and for a power-law exciting spectrum of the form 
J v = Jq^/^l)' 1 as 

x> _ Jo C dvfrA^V _ J„C(1 + x)kT 
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where ( denotes the integral. In terms of the ionization parameter S = P/J (, eq. (2.9) becomes 
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the transport equation (2.5) can also be simplified to give 
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where a VL = <j u (vl h ). Substituting for J and dJo/dz from eq. (2.9)-(2.11) we obtain the equation 
for the spatial behavior of the ionization fraction: 



where xq = x(z = 0) can be directly derived from eq.(2.10) 



Since both the integrand and the integral limits of eq. (2.12) depend on the temperature, some 
additional approximation is required. To estimate the value of the temperature we run the numer- 
ical code for an optically thin case to calculate the (constant) value of the temperature Tq to be 
substituted in eq. (2.13) as a function of E. In addition, in order to integrate eq. (2.12), we assume 
that the temperature inside the cloud remains approximately constant at the value Tq. Since most 
of the cooling mechanisms depend on the electron density n e , this simplification is fairly good for 



a relatively high (x £ 0.70) ionization fraction. In this case we find that the equation of state of 
the gas closely resembles a polytropic law 



Pocp (n/n+1) , n = -6.25, (2.14); 

the properties of polytropes with negative indexes have been studied by Viala & Horedt (1974). 
The solution x(z) of eq. (2.13) is then implicitly given by 

where C\ = logxo/l^o — 1| — (xq — 1) _1 and now x is calculated for T = T obtained numerically. 
Fig. 1 shows various solutions of eq.(2.15) for different values of the ionization parameter, for a 
spectral index 7 = 1.4. 

The main feature of the solutions is that, as one may expect, i c tends to become smaller for 
large values of E but the front is smoother. Also, for low values of S, the ionized zone can be as 
deep as several Kpc, and that means that in low pressure regions HVCs are not expected to have a 
dense, cold core. Although only approximate, the solutions found can be used to gain a qualitative 
insight of the problem and may serve as a test for the numerical scheme. 

In the following we will discuss the two main input parameters of the problem, namely the 
radiation field and the pressure of the Galactic halo. 

2. 1 Halo radiation field 

The ionizing radiation field in the halo is the sum of different components: the stellar radiation 
field from the disk, the quasar component and the X-ray background. As far as the disk contribution 
is concerned we will neglect it for the following reason. According to Brown & Gould (1970) (there 
is no appreciable difference at low energies with respect to the improved photoelectric cross-sections 
given by Morrison & Mc Cammon 1983) the energy at which half of the radiation is absorbed by 
gas in the Galactic disk is 

U; 1 (2.16) 



keV / \5.9 x 10 21 cm- 2 

If we consider the value of Nh = 3.1 x 10 20 cm -2 appropriate for the thick HI exponential distribu- 
tion as given by Dickey & Lockman (1990), then E* = 0.33 keV. Thus, unless the HI is very clumpy 
(a claim that has no actual observational support, Jahoda et al. 1985), it seems very unlikely that 
any source of hard photons located in the stellar disk may contribute sensibly to the halo radiation 
field. Photons from a population of hot stars like the central stars of planetary nebulae and white 
dwarfs with a larger scale height could be able to penetrate the layer, but their number is quite 
uncertain. In addition, even if present in appreciable number, these photons are barely sufficient to 
maintain the ionization of the electron component that is known to extend out to about ~ 1.0 kpc 
as pointed out by Reynolds (1990), (1993) and Nordgren, Cordes & Terzian (1993). In this sense 
the electron component almost totally insulate the halo from the Galactic ionizing radiation. The 
same conclusion can be drawn if one considers photons escaping from the open chimneys (Norman 
& Ikeuchi, 1989): again the majority of the photons are needed to maintain the ionization of the 
electron layer (Norman & Panagia 1991). Thus it appears that works including the contribution 
from photons escaping in the halo from the interior of multi-supernova remnants (Bregman & 
Harrington 1986) are overestimating the ionizing flux. 

A model for the radiation field in the halo proposed by Fransson & Chevalier (1985) has been 
widely used. The spectrum given by those authors, is the sum of the quasar component obtained by 



Sargent et al. (1979) using the Schmidt luminosity distribution and the extragalactic X-ray back- 
ground measured by Schwartz (1979). Lately, new studies have become available either theoretical 
and observational. While the X-ray portion of the spectrum given by Schwartz (1979) has been 
confirmed by recent observations by Ginga (for a review see Mc Cammon & Sanders 1990), for what 
concerns the QSO radiation field the newly obtained quasar luminosity functions have produced 
different estimates of the QSO contribution to the local UV extragalactic background (Terasawa 
1992, Madau 1992). In particular, Madau (1992) includes the effects of the opacity associated with 
the intervening Lya clouds and and Ly-limit absorption systems. Therefore, when these results are 
considered for the UV region of the EBR, we obtain the adopted spectrum 



5.96 x lO" 24 ^ 4 13.6ey < E < l.bkeV; 

5.1 x 10- 26 ^ e ° y 4 e- B/41 keV 1.5keV < E. (2.17) 



Note that this value is in agreement with the estimate I v {Ry) ~ 6 x 10 24 derived by Kulkarni & 
Fall (1993) from a proximity effect study in the distribution of Lya forest lines at low redshift. 

To take into account the possible radiative contribution from the hot confining gas, we have also 
used a composite spectrum which is a sum of the EBR and the hot gas emission contributions. For 
a plasma temperature 7\ ~ 10 6 K the continuum emission is dominated by the free-free processes 
(Landini & Monsignori Fossi 1990). The intensity of the free- free emission can be written as 

lfJ{E) = 1.66 x 10- 2 \EM)Tr°- 5 g(v,T h )Y ^e~ E / kT \ (2.18) 

n H 

where EM is the emission measure, Th the temperature of the plasma, g(u, Th) is the Gaunt factor, 
Z is the ion charge and n z /riH its abundance; following Landini & Monsignori Fossi (1990), we have 
included only H and He. The adopted values for the hot gas are taken from Burrows & Mendenhall 
(1991) who have derived T h = 1.25 x 10 6 K and EM = 0.006 cm" 6 pc. 

However, from an inspection of Fig. 9 of Landini & Monsignori Fossi (1990), it appears that, 
at the Th of interest, the total cooling losses are dominated by line emission, whose emissivity is 
about 18 times larger than the continuum one. In order to take into account this aspect, although 
in a rough manner, we assume that the spectral shape of the total emission from the hot gas 
(continuum+lines) is the same as the free-free but with a coefficient larger by a factor 18. This 
should be a reasonable assumption since the overwhelming majority of the lines for such a plasma 
are in the region shortwards of 912 A. 

2.2 Halo pressure 

The determination of the pressure the Galactic halo is a difficult task because it involves a firm 
understanding not only of the disk/halo interaction, but also of the processes occurring at the 
interface with the intergalactic medium. A great deal of uncertainty, as already discussed in the 
Introduction, persist on the nature (if any) of the disk/halo gas circulation. This topic has been 
recently reviewed by Spitzer (1990), who has also pointed out that the relationship among quiet 
phenomena as the buoyant gas in a fountain and highly dynamic transient ones is far to be clear. In 
addition, the question if the halo can be understood on average as a steady, and perhaps even static, 
structure still remains unanswered; nevertheless, many authors have investigated one-dimensional 
hydrostatic models for the halo (Fransson & Chevalier 1985; Hartquist & Morfill 1986; Bloemen 
1987; Boulares & Cox 1990). The purpose of this § is to use some of these results together with 
simple estimations to obtain reasonable limits to the halo pressure. 

2.2.1 Dynamical models 



Recently, Li & Ikeuchi (1992), have investigated the formation of giant halos around spiral 
galaxies. Their paper provides useful information about the pressure distribution in dynamically, 
fountain-dominated halos. Among the three types of halos they find (wind, bounded and cooled) 
only the last one (cooled) forms HVCs by thermal instabilites in a region between < w < 20 kpc 
and < z < 10 kpc, well out of the main galactic disk. The density and temperature vary quite 
rapidly in this region as can be realized from their Figs. 8(a)-8(b) and, therefore, rather than 
looking for a regular pattern we just derive an upper and lower limit to the pressure: 

1.3 x 10" 17 < P < 2.3 x 10" 14 ergs cm -3 ; (2.19) 

all the intermediate values are present in the cloud forming layer. This halo corresponds to a a 
temperature and density at the disk T = 10 6 K and no = 5 x 10~ 3 cm -3 , therefore log Po/k = 3.7 
It has to be pointed out that the pressure derived by Li & Ikeuchi is generally much lower than 
the one of the intergalactic medium (see next section), whose presence they did not consider. This 
could represent an inconsistency for their model. 

2.2.2 Static models 

The static model can be ideally divided in "hot" and "cold" ones, i.e., in which the gas support 
comes from its own thermal pressure or form pressure contained in other phases (magnetic fields, 
cosmic rays, turbulent motions), respectively. 

The simplest model for a hot halo is the isothermal one. The 2D pressure distribution in such 
a model can be derived from the hydrodynamic in a straightforward manner. The static Bernoulli 
equation can be written as 

(? a ]agp{m,z) + <t>{w,z)=il){w), (2.20) 

where <p is the Galactic gravitational potential, w is the galactocentric radius and ip is an integration 
constant with respect to z that can be obtained from a comparison with the tu-component of the 
steady momentum equation. For a constant galactic circular velocity v v , = v^\og{w) + C%, 

where Ci is an integration constant. Therefore 

P(w,z) =P e ( —Y e ^o-^,z))/c^ (2 21) 

Wo/ 

where (3 = v v /c s . In order to avoid an unphysical result for w — > oo we must assume that far 
away (3 — > 0. This leads to the final expression 

P(w, z) = p 0e (*0-*(~.*))/< (2.22) 

When the limit 

lim P(w, z) = P e* Q/c * = P IG 

z — >oo 

where Pig is the intergalactic pressure, is taken, given the values of Pig and P , c s is uniquely 
determined. The value of the gravitational potential at the solar radius is <I>0 = 4.8 x 10 14 cm 2 s -2 
for a circular velocity v v (wq) = 220 km s _1 (Binney & Tremaine 1987). Assuming for the in- 
tergalactic medium Tig = 10 8 K and mc = 10~ 6 cm -3 , and for P Q = 3.5 x 10~ 12 ergs cm -3 
(Cox 1990), the derived temperature is T ~ 10 6 K, exactly the same used as a boundary condition 
by the dynamical model of Li & Ikeuchi (1992). A plot of the pressure isocontour obtained from 
eq.(2.22) is shown in Fig. 2. The assumed potential is taken from de Boer (1990), also used by 
Ferrara (1993) to explain the vertical equilibrium of the Lockman component. It consists of a Oort 
vertical distribution and an exponential stellar radial potential: 



$(m,z) = <Jg{m) 



logcosh(z/z ) + ^e(w)(z/z ) 2 



(2.23) 



where z = 250 pc, a 2 g {w) = (15.4 km s" 1 ) 2 exp (w Q - w)/0A4w & , and e = 0.04,0.07,0.14 for 
w = 5, 10, 15 kpc, respectively. A quadratic interpolation to e has been used whenever required. In 
the region of interest (0 < w < 20 kpc, < z < 10 kpc) the lower and upper limits for a hot static 
model from Fig. 2 are 

10" 15 < P < 10" 12 ergs cm -3 (2.24). 

In the cold halo of models the support in the gravitational field is provided by sources different 
from thermal energy, as discussed above. We will refer, as a prototype, to the one-dimensional model 
presented by Boulares & Cox (1990), which appears to include the most updated data compilation 
of the distribution of the various phases of the ISM. A fit of the exponential component (which is 
the important one at high z) of the pressure from their results gives the following expression for 
the solar neighborhood 

P(z) = 2.61 x 10- 12 e- z/1 - 23kpc ergs cm" 3 . (2.25) 

We use this result in the same region considered for the Li & Ikeuchi model (3 < z < 10 kpc). In 
this region the lower and upper limits for the Boulares & Cox cold static model from eq.(2.25) are 

2.2 x 10 -13 < P < 7.7 x 10 _16 ergs cm" 3 . (2.26) 

The range of pressures obtained for the three different (dynamical, static hot/cold) types of 
halos is large and hence pretty unsatisfactory. This fact reflects the extremely poor comprehension 
of the halo and strongly recommends futher theoretical and observational studies. 

3. RESULTS AND THE DISTANCE METHOD 

The method we are proposing in order to determine the distance to the HVCs is essentially based 
on their detailed ionization structure, which is completely determined by the 5 coupled equations 
(2.1)-(2.5). For this purpose we have solved those equations numerically for each point inside the 
cloud, supposed to have a slab geometry and subjected to an external fixed pressure P. Figs. 
(3)- (5) show some of the results obtained for different values of P. The radiation field, illuminating 
the cloud from both sides, has been held fixed at this stage to the fiducial value given by eq. (2.17), 
and this particular choice will be referred as the standard field. The two panels of each figure show 
the density and temperature profiles as a function of the ratio z/£, where z is the depth inside the 
cloud and £ is the cloud size. Note that different curves refer to different values of the HI column 
density and, therefore, of £; obviously, solutions are symmetric about the point z/£ = given that 
the symmetry properties of the problem. 

As can be realized at a first glance to Figs. (3)-(5), the ionization structure of the cloud strongly 
depends on the pressure of the external medium. Following the discussion in § 2.2, we have selected 
three different values of the pressure P that presumably may well represent the actual values, 
namely P = 10 -15 , 10 -14 , 10 -13 ergs cm -3 . We have disregarded the value P = 10 -12 erg/cm 3 
because such a high pressure arises only in the static hot halo model and in regions close to the 
disk, where the contamination due to the galactic radiation field and local gaseous phenomena as 
outflows can become dominant. For P = 10~ 13 ergs cm -3 , the cloud is initially almost optically 
thin (£ = 0.6 pc, Nh = 2 x 10 16 cm -2 ) with an ionization fraction x ~ 0.42 and a temperature 
T ~ 9000 K. When the size is increased, self-shielding effects become more and more important, 
especially in the central parts of the cloud, and the ionization and the temperature decrease steadily 
due to the decreased number of ionizing photons available. The temperature structure will develop 
a maximum not located at the cloud edge, but slightly displaced in the interior; this is due to the 
fact that close to the edge, the temperature is determined by the larger availability of electrons for 



cooling rather than by a decrease in the number of ionizing photons. Eventually (last two curves), 
the temperature starts to fall down in an abrupt manner, while the ionization preserves its regular 
pattern. Continuing to increase Nh, the cloud develops a cold core, due to photon exhaustion in 
the partially ionized zone. This fact confirms our simple analytical estimates of § 2, shown in Fig. 
1. The characteristics of this central dense (n ~ 40 cm -3 ), cold (T ~ 20 K) core are not directly 
relevant to the aim of this paper (anyway, some discussion is given below). The important point 
is, however, the existence of a critical column density N H , corresponding approximately to the 
curves with the lowest ionization fraction shown in each figure, after which a central condensation 
completely shielded from the external radiation field is formed; in this case N H ~7x 10 17 cm -2 . 
We will exploit this fact, as explained below, to determine the distance to the HVCs. 

The same kind of behavior shown by the solutions for P = 10 -13 ergs cm -3 , can be qualitatively 
recognized for the two other pressure values adopted. For P = 10~ 14 ergs cm -3 the ionization frac- 
tion is higher, particularly when the medium is still optically thin [x ~ 0.8), and N H ~ 10 19 cm -2 . 
Finally, when the pressure is set to P = 10 -15 ergs cm -3 , the medium is almost completely ion- 
ized and in order to reach the critical column density the cloud size must be implausibly large 
(£ ^ 30 kpc) if compared with any realistic model of the HVCs. 

The previous results can be summarized in a Nh — I plane as the one presented in Fig. (6). 
Each series of points represents the Nh — £ relation for the curves previously discussed. According 
to the previous discussion, the plane can be divided in two subregions (Fig. 6): a lower one, in 
which the ionization front is "matter bounded" and clouds do not have any cold core, and an upper 
one in which the front is "radiation bounded" and a central cold core is found. Clouds located in 
the upper part of the plane have thus a composite structure made up by a neutral core surrounded 
by a warm envelope which constitutes an interface with the external, presumably hot, medium. A 
simple approximate expression for the "critical" curve N H {£) is found to satisfy 



with Ni = 3 x 10 17 cm 2 . An analogous relation, which is a very good approximation in the range 
5 x 10 15 ^ P ^ 10~ 13 ergs cm -3 , between N H and P is 



where N 2 = 7.4 x 10 17 and P-13 = P/10 -13 . In addition to the case in which only the EBR field 
is considered, the analougous results for the sum of the EBR + free-free ionizing spectrum is also 
reported in Fig. 6. This addition is not modifying the results in a substantial way because of the 
relatively low EM of the hot gas, and hence of the low intensity of its radiation field. Note that there 
are no open triangles in the case P = 10 -15 ergs cm -3 because they are physically meaningless, 
indicating cloud sizes of the order of 10 2 kpc, larger than the halo itself; one may wonder about the 
extragalactic nature of HVCs (which is excluded by almost any sort of data), but the hypothesis 
of this paper should then be radically modified. 

It is natural to identify the core+interface structure found with the observational evidences 
reported by many authors and discussed in the Introduction, even if not much significance should 
be attached to the precise values of the line widths since non thermal mechanisms may concur to 
the line broadening; the important point is the observed values, typically ^ 20 km s _1 , are higher 
than the ones corresponding to our solutions. At this point the suggested strategy to determine 
the distance to the HVCs becomes almost self-explanatory. In order to be as clear as possible we 
put it in a schematic form: 

1.- Using a cloud with a core+interface structure, measure its Nh in the interface; the measured 
Nh will hence correspond to the critical column density N H ; 




(3-1) 



N H (£) = N 2 P:^ cm 



-2 



(3.2) 



2. - The value of N H must be located on the critical curve (3.1). This gives the linear size of the 

interface, which, given the existence of the core, must be exactly the £ corresponding to that 
critical column density. By the way, pressure can be therefore immediately determined through 
relation (3.2); 

3. - Once the linear size of the interface £ is known, from its angular diameter 9, the distance A can 

be found through the relation A = £/29. 

The procedure is quite straightforward if the cloud is observed to have the required core+interface 
structure and column density measures are available. Note that the actual characteristics and 
column density of the core are not relevant to the distance determination. The mere existence 
of the core bounds the column density of the partially ionized zone (interface) to the appropriate 
critical density for the given pressure. Some inferences on the physical state of the cores can 
nevertheless be made. The very low temperature characterizing the cores in our model is to a large 
extent determined by the only heating mechanism included apart from the radiative one, i.e., C 
ionization. Other heating mechanisms can of course be foreseen; however, the non-detection of 
any HVCs in the IRAS maps (Wakker & Boulanger 1986), if interpreted as a low temperature of 
the dust, favours the idea that HVCs are a cool environment. However, the non-detection can be 
alternatively explained as due to a low dust abundance, since dust can be heavily sputtered during 
the expulsion process from the disk (Ferrara et al. 1991). 

Adding another step to our procedure, we can also infer an approximate temperature of the core: 

4. - Using the measured core Nh, from the determination of its linear size £ c which is possible 

knowing A, we obtain an average value of the density in the core (n c ), from which T c can be 
derived from the equation of state. 

From a theoretical point of view, there are two main mechanisms which can introduce some error 
in the distance determined according to the scheme presented. The first one is thermal conduction 
(Cowie & McKee 1977; Draine & Giuliani 1984; Begelman & Mc Kee 1990; McKee & Begelman 
1990; Ferrara & Shchekinov 1993) whose main effect is to broaden the ionized zone. However, 
if clouds are formed from a thermal instability in the fountain, by definition their size must be 
larger that the critical wavelength at which thermal instabilities become stabilized by conduction 
(Field 1965). Ferrara k, Shchekinov (1993) have demonstrated that the dynamical effects of thermal 
conduction in such conditions are negligble. Also conductive interfaces at the steady state tend to 
have a remarkably flat temperature profile: for a sperical cloud, T(r) = Th(l — £/r) 2 ^ 5 , where r is 
the distance from the cloud center and r > £, and thus most of the interface is at a temperature close 
to the one of the hot medium T^. In these condition virtually all the hydrogen is completely ionized 
and the contribution of this gas to Nh is negligible. A second mechanism providing ionization is 
mechanical input by shocks. The sound speed in the hot gas is 



where it has been assumed that the mean molecular weight fi = 0.65, m p is the proton mass, and 
n_3 = ra/j/10 -3 is the hot gas density; thus, by definition, some HVCs move subsonically for the 
reference values of eq. (3.3). The emission measure integrated on a path of 10 kpc for such a hot 
gas is E m (ri-3 = 1) = 10~ 2 cm -6 pc, consistent with the one derived by the ROSAT soft x-ray 
shadowing experiments (Burrows Sz Mendenhall 1991; Snowden et al. 1991); also the temperature 
derived by a fit to the data is consistent with the estimate (3.3). Therefore, it appears that a 
relevant fraction of the HVCs should not be affected by a shock, unless a large transversal velocity 
component is present. However, shocks may have some importance in the thermal balance of the 
core, and eventually may lead to episodes of star formation, according to the calculations made by 




(3-3) 



Dyson & Hartquist (1983). Viscous heating, TL oc phV 3 is not going to be important given the low 
density of the hot medium. We conclude that the error on the distance determination due to the 
effect of thermal conduction is negligible; in addition, for those HVCs with velocity lower than c s 
specified by eq. (3.3), shocks are not expected to appear. 

Probably more important are the possible errors related to real data. It could be difficult to 
locate observationally the position of the inner and outer edges of the interface in an accurate 
manner. The inner edge is defined as the limit of the region with narrow line profiles; however, 
this edge is often not regular. For the outer edge, instead, it is crucial to resort to high sensitivity 
data to detect the low HI column density associated with the external regions of the interface. In 
addition, our assumption of a slab geometry can introduce some additional error, particularly when 
the interface has dimensions much larger than the core. This could be avoided by appropriately 
selecting the cloud sample. With respect to the last point, we mention that in the framework of the 
fountain model, HVCs are predicted to have a a sheet-like shape (Kahn 1991) which would make 
the slab geometry very appropriate. 

From the previous discussion it appears that the direct application of the proposed method 
requires particular attention to both sample selection and the estimation of possible errors. We 
plan to present these results in a forthcoming publication; we stress, however, that a dedicated 
observational study could serve this purpose much better. 



In the previous section we have presented a simple method based on the photoionization structure 
of the HVCs which can be used in order to obtain reliable distances to the clouds possesing a two- 
phase structure composed by a warm ionized interface and a cold neutral core. However, our 
treatment provides an opportunity to have some insight on some global properties of the Galactic 
halo and the EBR. This is posssible if we compare the Ha emission as predicted by our models 
with the one detected by Kutyrev & Reynolds (1989) (hereafter KR) and Songaila, Bryant & Cowie 
(1989) (hereafter SBC). The first authors report the detection of Ha emission from a HVC situated 
in Cetus with a surface brightness I a = 8.1 ± 1.9 x 10 -2 R (where one rayleigh corresponds to 
10 6 /47r photons cm -2 s _1 sr _1 ); Songaila et al. deduce instead I a = 3 x 10~ 2 R from two points in 
a HVC belonging to the Complex C. These are the only detections we are aware of and therefore, 
given the restricted sample, no definite conclusions can be drawn; the comparison could be anyway 
of some interest. 

Assuming Case B recombination, the Ha intensity integrated on our model slab cloud can be 
written as 



where K = 8.17 x 10 -16 , n is the total gas density, x is the fractional ionization and h c ,IiHe ar e 
the abundances of C and He relative to H, respectively. The results obtained applying eq. (4.1) to 
the three models for low, intermediate and high pressure cases discussed in the previous section are 
shown in Fig. 7 along with the oservational results. It is immediately realized that the observed Ha 
intensity is larger than the theoretical one for all the pressure cases considered and the discrepancy 
becomes larger as the pressure of the external medium is increased. The explanations for this 
dicotomy can all be found in two basic possibilities, namely i) the intensity of the radiation field 
must be larger than the one assumed, and/or ii) some mechanical energy input must take place. As 
for the first point, it is difficult to conceive a dramatically different value from the EBR intensity 
given by eq. (2.17): in fact, it has been proven able to explain several effects at high redshift and 
confirmed by observations in its high energy part. If we exclude that any ionizing photon from 



4. FURTHER IMPLICATIONS 




(4.1) 



the galaxy can escape the thick Reynolds layer, the only remaining source of ionizing photons may 
come from free-free emission of the hot gas confining the cloud. In the previous Section we have 
seen that some modification to the ionization structure of the cloud and to the interface size can 
be introduced by the consideration of this contribution. However, the Ha emission calculated from 
the models including the radiation field from the local hot gas show an enhancement only of a few 
percent with respect to the pure EBR field case and certainly not sufficient to explain the observed 
excess. 

Thus we are left with the second possibility, i.e., Ha emission is provided by recombinations 
following a shock wave created by the HVC/hot gas interaction. Support for this hypothesis comes 
from the two following evidences. First, the two clouds have velocities well in excess of the sound 
speed in the hot medium (eq. [3.3]): for the KR cloud \vlrs\ ~ 300 km s _1 , while for the SBC 
one IvLiJsl ~ 140 km s -1 ; second, the Ha emission excess is correlated with the velocity as one 
would expect if part of the ionization is produced by a shock. We conclude that the suggested 
method is not applicable to those clouds in a straightforward manner since we have assumed that 
the ionization of the clouds is due to the radiation field (EBR+free-free) only; at any rate the 
distance method in such cases will provide at least a lower limit to the distance. This could be 
particularly useful if absorption features towards a background star are detected from the cloud. 
Finally, the calculated Ha emission (Fig. 7) provides a powerful test of the conditions in which the 
proposed distance method can give reliable results, and clearly indicates the possible presence of a 
shock. 
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FIGURE CAPTIONS 



Figure 1 Ionization structure of the cloud partially ionized zone. The curves show the solutions 
of eq. (2.15); the numbers refer to different values of the ionization parameter S = P/Jo(; z is 
measured from the cloud edge. 

Figure 2 Pressure isocontours for a static, isothermal, hot halo; the temperature and pressure at 
the disk are T = 10 6 K, P Q = 3.5 x 10" 12 erg cm" 3 . 

Figure 3 Ionization and temperature structure for a cloud of size I exposed to the standard EBR 
field given by eq. (2.17); the pressure is constant at a value P = 10~ 13 erg cm -3 . Each line refers to 
a different value of the column density Nh which is in the range 3 x 10 16 < Nh < 7.5 x 10 17 cm -2 . 

Figure 4 Same as Fig. 3, but for P = 10" 14 erg cm" 3 , and 5 x 10 16 < N H < 8.5 x 10 18 cm" 2 . 

Figure 5 Same as Fig. 3, but for P = 1CT 15 erg cm" 3 , and 1 x 10 15 < N H < 1.4 x 10 18 cm" 2 . 

Figure 6 Hydrogen column density of the cloud as a function of its size; the numbers refer to 
different values of the pressure in erg cm -3 . Solid triangles refer to the standard EBR field given 
by eq. (2.17), open triangles are for the composite field EBR+free-free from hot gas. The upper 
dashed part describes the region of the parameter space in which the clouds develop a core+interface 
structure in the EBR field. 

Figure 7 Calculated Ha intensity, I a , from the cloud interface as a function of cloud size for the 
three models of Fig. 6; the numbers refer to different values of the pressure in erg cm -3 . Also 
shown (horizontal lines) are the observed values for two HVCs as given by RK and SBC. 



